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The DASI discovery of CMB polarization has opened a new chapter in cosmology. Most of the 
useful information about inflationary gravitational waves and reionization is on large angular scales 
where Galactic foreground contamination is the worst, so a key challenge is to model, quantify and 
remove polarized foregrounds. We use the POLAR experiment, COBE/DMR and radio surveys to 
provide the strongest limits to date on the TE cross power spectrum of the CMB on large angular 
scales and to quantify the polarized synchrotron radiation, which is likely to be the most challenging 
polarized contaminant for the MAP satellite. We find that the synchrotron E- and B-contributions 
are equal to within 10% from 408 — 820MHz with a hint of i7-domination at higher frequencies. We 
quantify Faraday Rotation & Depolarization effects in the two-dimensional {i, i/)-plane and show 
that they cause the synchrotron polarization percentage to drop both towards lower frequencies and 
towards lower multipoles. 



I. INTRODUCTION 

The recent discovery of cosmic microwave background 
(CMB) polarization by the DASI experiment [1] has 
opened a new chapter in cosmology - see Figure 1. Al- 
though CMB polarization on degree scales and below can 
sharpen cosmological constraints and provide important 
cross-checks [2,3], the potential for the most dramatic 
improvements lies on the largest angular scales where it 
provides a unique probe of the reionization epoch and pri- 
mordial gravitational waves. For instance, forecasts [4,5] 
indicate that the MAP satellite can measure the reion- 
ization optical depth r seventeen times more accurately 
using polarization information, and that polarization in- 
creases the sensitivity of the Planck satellite to tensor 
modes by a factor of 25. 

Unfortunately, these large scales are also the ones 
where polarized foreground contamination is likely to be 
most severe, both because of the red power spectra of dif- 
fuse Galactic synchrotron and dust emission and because 
they require using a large fraction of the sky, including 
less clean patches. The key challenge in the CMB pola- 
rization endeavor will therefore be modeling, quantifying 
and removing large-scale polarized Galactic foregrounds. 
This is the topic of the present paper. We will use the 
POLAR experiment to provide the strongest limits to 
date on cross-polarized microwave background and fore- 
ground fluctuations on large angular scales, and employ 
polarization sensitive radio surveys to further quantify 
the polarized synchrotron radiation, which is likely to 
be the most challenging contaminant in the polarization 
maps expected from the MAP satellite. 
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FIG. 1. Summary of constraints on polarization so far. From 
top to bottom, the three curves show the concordance model predic- 
tions for Cj, Cf and , respectively. Four reionization models 
with r=0.1, 0.2, 0.3 and 0.4 are also plotted (left dotted lines from 
bottom to top in both plots). The limits for E are shown in the 
upper panel: Penzias & Wilson 65 [6], Caderni 78 [7], Nanos 79 
[8], Lubin & Smoot 79 [9], Lubin & Smoot 81 [10], Sironi 98 [11], 
Lubin 83 [12], SASK (W93 [13],N97 [14]), TOCO (T99 hexagons 
[15]), P88 [16], F93 [17], P97 [18], SCO [19], DMR [20], PIQUE 
(H02 [21]) and POLAR (KOI [22]). The limits for X are shown in 
the lower panel: PIQUE (dOC02 [23]) and POLAR ("This Work"). 
The shaded regions are the DASI results (Kv02 [1]). 
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At microwave frequencies, three physical mechanisms 
are known to cause foreground contamination: syn- 
chrotron, free-free and dust emission. When coming from 
extragalactic objects, this radiation is usuaUy referred to 
as point source contamination and affects mainly small 
angular scales. When coming from the Milky Way, this 
diffuse Galactic emission fluctuates mainly on the large 
angular scales that are the focus of this paper. Except for 
free- free emission, all the above mechanisms are known to 
emit polarized radiation. In the near term, the best mea- 
surement of large-scale polarization will probably come 
from the MAP satelhte. At MAP's frequency range (22- 
90 GHz), synchrotron radiation is likely to be the dom- 
inant polarized foreground [4]. Unfortunately, we still 
know basically nothing about the polarized contribution 
of the Galactic synchrotron component at CMB frequen- 
cies [4,24-29], since it has only been measured at lower 
frequencies and extrapolation is complicated by Faraday 
Rotation. This is in stark contrast to the CMB itself, 
where the expected polarized power spectra and their de- 
pendence on cosmological parameters has been computed 
from first principles to high accuracy [30-33] . 

Polarization of the Galactic continuum emission was 
first clearly detected in 1962 [34]. In the succeeding 
years, polarization measurements of the northern sky 
were made at frequencies between 240 and 1415 MHz 
(see [35] and references therein) with resolutions of only 
a few degrees. No large-area survey has been published 
since the compendium of Brouw and Spoelstra [36] and 
high-resolution surveys have only begun to be made re- 
cently. The first major investigation done after [36] is 
that of [37] , who observed a section of the Galactic plane 
defined by 49° < ^ < 76° and \b\ < 15°, at frequency of 
2.7 GHz. The study of [38] provides the highest resolu- 
tion insight into the small-scale structure of the Galaxy; 
however, this only covered a few areas of the sky which 
were not larger than a degree or so across. Recently, two 
fully-sampled polarimetric surveys were done at 2.4 GHz 
[39,40] and 1.4 GHz [41,42]. AU of these high-resolution 
surveys covered only regions near the Galactic plane, so 
in order to use them for inferences relevant to CMB ex- 
periments, they need to be extrapolated both in Galactic 
latitude and in frequency. 

The rest of this paper is organized as follows. In section 
Section II, we review the basics of CMB and synchrotron 
polarization as well as our methods for measuring and 
modeling it. We present our results in Section HI and 
discuss our conclusions in Section IV. 
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FIG. 2. Examples of CMB polarization, showing how the reion- 
ization optical depth t affects the T and E power spectra (top) 
and the TE correlation ri (bottom). Solid, dashed and dotted 
curves correspond to t=0, 0.2 and 0.4, respectively. As discussed 
in [23], changing the cosmological parameters affects the polarized 
and unpolarized power spectra rather similarly except for the cases 
of reionization and gravitational waves. In the reionization case, 
a new series of peaks are generated at large scales. Top-panel: 
Although there is no visible change in T at large scales, there is 
clearly a visible change in E since the Sachs- Wolfe nuisance is un- 
polarized and absent. Lower-panel: On small scales, reionization 
leaves the correlation unchanged since Cj" and are merely 
rescaled. On very large scales, ri drops since the new polarized 
signal is uncorrelated with the old unpolarized Sachs-Wolfe signal. 
On intermediate scales £ > 20, oscillatory correlation behavior is 
revealed for the new peaks. For more details about CMB polariza- 
tion and reionization see [43] 



II. PHENOMENOLOGY 
A. Notation 

CMB measurements can be decomposed into three 
maps {T,E,B), where T denotes the unpolarized and 
{E,B) denote the polarized components, respectively. 
Note that an experiment that is insensitive to polariza- 
tion does not measure T but rather that total (unpolar- 
ized plus polarized) intensity; although this distinction 
has traditionally been neglected for CMB experiments 
where the polarization fraction is small, it is important 
both for foregrounds (which can be highly polarized) and 
for precision CMB experiments. From these three maps 
we can measure a total of six angular power spectra, here 
denoted by Cj , CF , Cf , Cf , CY and , correspond- 



2 



ing to the TT, EE, BB, TE, TB and EB correlations*, 

respectively. 

By parity, Cf — Cf = for scalar CMB fluctuations, 
but it is nonetheless worthwhile to measure these power 
spectra as probes of both exotic physics [44-46] and fore- 
ground contamination. Cf = for scalar CMB fluctua- 
tions to first order in perturbation theory [30-33] - sec- 
ondary effects such as gravitational lensing can create B 
polarization even if there are only density perturbations 
present [47] . In the absence of reionization, Cj is typi- 
cally a couple of orders of magnitude below Cf on small 
scales and approaches zero on the very largest scales. 

The cross-power spectrum C-^ is not well suited for 
the usual logarithmic power spectrum plot, since it is 
negative for about half of all values [23] . A theoretically 
more convenient quantity is the dimensionless correlation 
coefficient 

„X - W /IN 

^ " (CjCf)V2' 

plotted on a linear scale in Figure 2 (lower panel), since 
the Schwarz inequality restricts it to lie in the range 
— 1 < rf < 1^ . Prom here on we use re as shorthand 
for rf . For more details about and how it depends on 
cosmological parameters, see section II. b in [23]. 

B. Our Knowledge of Synchrotron Emission 

The Galactic InterStellar Medium (ISM) is a highly 
complex medium with many different constituents inter- 
acting through a multitude of physical processes. Free 
electrons spiraling around the Galactic magnetic field 
lines emit synchrotron radiation [48], which can be up 
to 70% linearly polarized (see [49,50] for a review). 

The power spectrum Ce of synchrotron radiation is 
normally modeled as a power law in both multipole £ 
and frequency u, which we will parametrize as 

ST^{u)=A(-j with A^u^", (2) 

where STe = [£{£ + l)Ce/2Tr]^/^. This definition implies 
that Ce oc i''' for ^ 3> 1 and that the fluctuation ampli- 
tude oc . The standard assumption is that the total in- 



*Prom here on, we adopt the notation TT = T, EE = E, 
BB = B,TE = X,TB = Y and EB = Z. 

^Note that for experiments where CMB polarization is mea- 
sured with a very low signal-to-noise ratio, Cf is a more useful 
quantity than rf . This is because they may be able to place 
upper and lower limits on Cf but can place no meaningful 
limits on rf unless they can statistically rule out that Cf in 
the denominator of equation (1). 



tensity has a « —2.8 with variations of order 0.15 across 
the skyt [51]. 

As to the power spectrum slope /?, the 408 MHz 
Haslam map [55,56] suggests /3 of order -2.5 to -3.0 down 
to its resolution limit of ~ [58-61]. A similar analysis 
done on the 2.3 GHz Rhodes map of resolution 20' [53] 
gives P = -2.92 ± 0.07 [62] (flattening to /3 « -2.4 at 
low Galactic latitudes [29]). 

For the polarized synchrotron component, our obser- 
vational knowledge is, unfortunately, not as complete. 
To date, there are measurements of the polarized syn- 
chrotron power spectrum obtained basically from three 
different surveys [63]: the Leiden surveys** [36,35], the 
Parkes 2.4 GHz Survey of the Southern Galactic Planed 
[39,40], and the Medium Galactic Latitude Survey-'--'- 
[41,42,64]. 

These measurements exhibit a much bluer power spec- 
trum in polarization than in intensity, with /3 in the range 
from 1.4 to 1.8 [4,24-29]. These results are usually taken 
with a grain of salt when it comes to their implications 
for CMB foreground contamination, for three reasons: 

1. Extrapolations are done from low to high galactic 
latitudes; 

2. Extrapolations are done from low to high frequen- 
cies; and 

3. Much of the available data is undersampled. 

The Leiden surveys extend to high Galactic latitudes and 
up to 1.4 GHz but are unfortunately undersampled, while 



'Because the spectral index a depends on the energy dis- 
tribution of relativistic electrons [48], it may vary somewhat 
across the sky. One also expects a spectral steepening to- 
wards higher frequencies, corresponding to a softer electron 

spectrum ( [52]; Fig 5.3 in [53]). A recent analysis done at 
22 MHz [54] shows that a varies slightly over a large frequency 
range. 

^Although the interpretation is complicated by striping 
problems [57]. 

**The observations done by Brouw and Spoclstra covered al- 
most 40% of the sky extending to high Galactic latitudes. Us- 
ing the same instrument, they observed the polarized Galaxy 
in Q and U in five frequencies from 408 MHz up to 1.4 GHz 
and with angular resolutions from 2.3° at 408 MHz up to 
0.6° at 1.4GHz. Unfortunately this data was also undersam- 
pled, making it difficult to draw inferences about its polarized 
power spectrum. 

^^This survey covers a strip 127° long and at least 
10° wide centered in the Galactic plane, with a res- 
olution of FWHM=10.4'. It is pubUcally available at 
http:/ / www. uq. edu. au/^ roy/. 

**The Medium Galactic Latitude Survey maps the Galac- 
tic plane within ±20°, with a resolution of FWHM=9.35' 
at 2.4 GHz. This survey is partially available at 
http. -//www. mpifr-bonn. mpg. de/ staff /huyaniker /index, htm. 
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the Parkes and the Medium Galactic Latitude Surveys 
only probe regions around the Galactic plane. In the 
following three sections, we will discuss these three prob- 
lems in turn. 



1. The Latitude Extrapolation Problem 

Although only high Galactic latitudes are relevant for 
CMB work, most of the data used for understanding 
the polarized CMB foreground contamination is at low 
Galactic latitudes. Figure 3 shows that whereas the total 
intensity of the synchrotron emission depends strongly on 
the Galactic latitude, the polarized component is approx- 
imately independent of Galactic latitude — indeed, in the 
three polarized images, it is difficult to distinguish the 
galactic plane at all. As noticed long ago by [40] , there is 
a faint, quasi- uniform polarized component of the Galac- 
tic polarized emission in their survey, upon which the 
emission from other features is superimposed: towards 
the higher latitudes, this faint component appears simi- 
lar in both structure and intensity to the correspondent 
lower latitude emission. This well-know empirical result 
can be also seen (in a more quantitative way) in the Lei- 
den surveys. Figure 4 shows that in the frequency range 
between 408 MHz to 1.4 GHz, the Polarization Intensity 
P {P=^J + U'^) is basically constant as the Galactic 
latitude |6| increases, whereas the polarization insensitive 
surveys (such as the 408 MHz Haslam and the 1420 MHz 
Reich & Reich [65]) have the bulk of their emission com- 
ing from the Galactic plane. 

The usual interpretation of this very weak latitude de- 
pendence of polarized synchrotron radiation is that the 
signal is dominated by sources that are nearby compared 
to the scale height of the Galactic disk, with more distant 
sources being washed out by Depolarization (to which 
we return in the next subsection). As a result, having 
well-sampled polarized maps off the galactic plane at the 
same frequencies would not be expected to affect our re- 
sults much, since they would be similar to those in the 
plane. This issue, however, deserves more work as far 
as extrapolation to CMB frequencies is concerned: the 
latitude dependence may well return at higher frequen- 
cies as Depolarization becomes less important, thereby 
revealing structure from more distant parts of the Galac- 
tic plane. In this case, extrapolating from an observing 
region around the Galactic plane to higher latitudes may 
well result in less small-scale power in the angular distri- 
bution. 

If we are lucky, many of the complications of extrapo- 
lating to higher latitude may largely cancel out the com- 
plications of extrapolating to higher frequency, thereby 
making it easier to quantify the polarized CMB fore- 
ground problem. The reason for optimism is the follow- 
ing: at high latitudes (which is all that really matters 
for CMB research) , the foreground signal will be entirely 
due to nearby emission within the scale height of the 



thick Galactic disk; and at low frequencies in the Galac- 
tic plane (which is where we have really good data), the 
polarized signal we see may well be dominated by such 
nearby emission, with emission from more distant regions 
in the Galactic disk hidden by Depolarization. 




FIG. 3. The nature of the Galactic synchrotron emission. 
Clockwise from top left, the panels show Stokes T, (7, Q, and P 
(defined as P=^Q'2 + U'^) from Block 3 of the Parkes 2.4 GHz 
Survey of the Southern Galactic Plane. 
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FIG. 4. The polarized and total synchrotron component as a 
function of the Galactic latitude. Each of five Leiden polarized 
surveys was divided in six slices of eq ual area, we then calculated 
the mean intensity (defined as P=\/ + U^) for each of those 
slices. A similar procedure was used for the polarization insensi- 
tive Haslam and Reich & Reich surveys, but 12 slices were chosen 
instead. The top panel show the results from the five Leiden sur- 
veys plus the 408 MHz Haslam data, while the bottom panel show 
the results from the Leiden 1.4GHz survey and 1.42 GHz Reich &; 
Reich data. Comparison between polarized components and the 
total intensity at the same frequency illustrates that the polarized 
synchrotron is almost independent of the Galactic latitude while 
the unpolarized emission is strongly concentrated in the Galactic 
plane. 
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2. Faraday Rotation, Depolarization and the Frequency 
Extrapolation Problem 

The plane of a polarized wave may be regarded as the 
sum of two circularly polarized components of opposite 
handedness. In an ionized medium with a non-zero mag- 
netic field, these two components propagate with differ- 
ent phase velocities, which will result in a rotation of the 
plane of polarization of the linearly polarized radiation. 
This rotation, known as the Faraday Rotation^ , produces 
a change in polarization angle A9 of 

A6» = O.SIA^ / UeBiidL = X^RM [rad], (3) 

JO 

where A is the wavelength given in meters and the quan- 
tity M/X^ is called the Rotation Me asurc {RM — usually 
expressed in units of rad m~^). The integral is done over 
the line of sight from us to the emitting region at a dis- 
tance L in pc, He is the the free electron density in cm~^ 
and is the magnetic field parallel to the line of sight 
in /zG. 

From equation (3) it is easy to sec that observations of 
this synchrotron radiation in several frequencies allows 
the determination of Rotation Measures in the diffuse 
radiation. From the obtained structure in the Rotation 
Measure on different scales, we can obtain information on 
the magnetic field parallel to the line of sight, weighted 
with electron density - an example of this method can 
be found in [68]. In radio astronomy, Faraday Rotation 
has become one of the main tools to investigate the in- 
terstellar magnetic field (see, e.g. [69,70]). 

It is important to point out, however, that Faraday 
Rotation can only change the polarization angle and not 
the polarized intensity P. The fact that we do see struc- 
ture in P that is not correlated with a counterpart in 
intensity T implies that part of the radiation has been 
depolarized [38] . A simple visual comparison of the total 
intensity and polarized maps of the same region in the 
sky of the Parkes 2.4 GHz survey shows Depolarization 
at work (see Figure 3): many sources which present an 
intense total emission do not show a counterpart in the 
polarized maps; similarly bright regions of extended po- 
larization are not connected with unpolarized sources. A 
more detailed study of this same survey reached similar 
conclusions: Giardino et al. [29] showed that the E and 
B power spectra were dominated by changes in the po- 
larization angle rather than by changes in the polarized 
intensity, suggesting that Faraday Rotation was playing 
a significant role-l-. 



detailed discussion of the Faxaday Rotation & Depolar- 
ization effects as well as their importance in astrophysical 
observations is in given in [67]. 
* Although at first glance the images in Figure 3 suggest that 



Depending on the frequency and beamwidth used, De- 
polarization can play an important role in polarization 
studies of the Galactic radio emission [35]. As discussed 
by Cortiglioni and Spoelstra [71], Depolarization can 
have four causes: 1) differential polarization along the 
line of sight, 2) differential polarization across the beam, 

3) differential Faraday Rotation across the beam, and 

4) differential Faraday Rotation and polarization across 
the bandwidth. If the bandwidth is very narrow, wc can 
neglect item 4', also, if the polarized data has been suf- 
ficiently sampled, smoothing it to a largest beam may 
inform us about items 2 and 3, leaving us with item 1 as 
the expected main source of Depolarization^. 

Because of the complicated interplay of these mech- 
anisms, we should expect both the amplitude and the 
shape of the polarized synchrotron power spectrum to 
change with frequency. We will therefore take an empir- 
ical approach below and use the available data to map 
out (for the first time) the two-dimensional region in the 
(£, z/) plane where Faraday Rotation & Depolarization are 
important. 

3. Incomplete Sky Coverage and the Undersamplmg Problem 

For the case of undersampling in the Leiden surveys, 
some authors have overcome this problem by doing their 
Fourier analysis over selected patches in the sky where 
they believe the average grid space in the patch is close 
to the map's beam size, so that they can apply a Gaussian 
smoothing on it - this is well explained and illustrated in 
[27]. Fortunately, we can eliminate this problem by mea- 
suring the power spectra with the matrix-based quadratic 



the polarized and unpolarized components arc uncorrclatcd, 
[40] found that for some patches in their images there is a 
good correlation between the polarized and total power in- 
tensities. Therefore they conclude that a good fraction of the 
polarized omission seen over the plane was caused by changes 
in synchrotron cmissivity, rather them any Depolarization or 
Faraday Rotation of the synchrotron background. According 
to [40], variations in synchrotron emission can be caused by 
increases in the density of rclativistic electrons (due to SNRs), 
and/or variations in the magnetic field intensity. 

It is important to point out that the relative importance of 
these two mechanisms (Faraday Rotation & Depolarization 
and changes in the synchrotron emissivity of the source re- 
gions) over the Galactic plane region are currently unknown 
[66]. 

^In the case of Leiden surveys, item 4 is negligible. Based 
in previous analysis done over the Galactic loops at 1.4 GHz 

[72,73], Spoelstra [35] argued that items S and 5 have a rela- 
tively minor contribution to the Depolarization in those sur- 
veys. Leaving, therefore, differential polarization along the 
line of sight as the main source of Depolarization. 
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estimator technique that has recently been developed for 
analyzing CMB maps [74,75,23]. 

Although the undersampling and partial sky coverage 
results in unavoidable mixing between different angular 
scales £ and polarization types {E and B), this mixing 
(a.k.a. leakage) is fully quantified by the window func- 
tions that our method computes [75] and can therefore 
be included in the statistical analysis without approxi- 
mations. Specifically, wc compute the six power spectra 
{T,E,B,X,Y,Z) described in Section II A so that the 
leakage, if any, is minimal. 

In [75] it was argued that susceptibility to systematic 
errors could be reduced by choosing the "priors" that 
determine the quadratic estimator method to have van- 
ishing cross-polarizations, X = Y = Z = and it was 
shown that this simplification came at the price of a very 
small (percent level) increase in error bars. In Appendix 
A of [23], it was shown that this choice has an important 
added benefit: exploiting a parity symmetry, it elimi- 
nates 14 out of the 15 leakages, with only the much dis- 
cussed [75,32,76-79] E — B leakage remaining. In [80] it 
was shown that even the remaining E — B leakage can, 
in principle, be removed. Unfortunately, this technique 
cannot be applied here, since it works only for a fully 
sampled two-dimensional map. 



Table 1 - POLAR-DMR Power Spectrum 

Zcff zb M ST"^ ± a {[iK'^\ [/i-g]^"^ 

T 15.6±6.6 487.0±270.6 22.1+^-| 

E 12.6±4.5 -9.9± 32.0 <4.7 (7.4) 

B 12.6±4.5 13.9± 32.0 <6.8 (8.8) 

X 14.0±4.8 -26.0± 48.5 <8. 7(11.1) 

Y 14.0±4.8 -0.1± 48.5 <7.0 (9.8) 

Z 11.4±2.9 -50.0± 31.6 <6.6(10.7) 

'"'Values in parentheses are 2-a upper limits. Cross- 
correlation upper limits refer to jXj, jyj and \Z\. 



1. Quadratic Estimator Analysis 

We measure the six power spectra described in Sec- 
tion II A using the quadratic estimator method exactly 
as described in [75]. We computed fiducial power spec- 
tra with the CMBFAST software [84] using cosmological 
parameters from the concordance model from [85] (that 
of [86] is very similar). Table 1 shows the result of our 
band-power estimation. The values shown in parentheses 
in the rightmost column of this table are our 2-a upper 
limits. In these calculations, we used 5 multipole bands 
of width Ai = 6 for each of the six polarization types 
{T,E,B,X,Y,Z), thereby going out to £ = 30, and av- 
erage the measurements together with inverse-variance 
weighting into a single number for each polarization type 
to minimize noise. 

We used the combined DMR 53-^90 GHz data to ob- 
tain good sensitivity to the unpolarized component. We 
perform our analysis using strips of the DMR data of 
width ±15° around the POLAR declination - we found 
that further increasing in the width of these disks did not 
significantly tighten our constraints. Finally, we elimi- 
nated sensitivity to offsets by projecting out the mean 
(monopole) from the T, Q and U maps separately. 

The detection of unpolarized power is seen to be con- 
sistent with that published by the DMR [83] group. Ta- 
ble 1 shows that we detect no polarization or cross- 
polarization of any type, obtaining mere upper limits, 
just as the models predict. The window functions reveal 
substantial leakage between E and B, so that the lim- 
its effectively constrain the average of these two spectra 
rather than both separately. This large leakage is due 
to the one-dimensional nature of the POLAR dataset, 
and can be completely eliminated with a fully sampled 
two-dimensional map [80]. 

Finally, we perform the same analysis described above 
by replacing the DMR stripe with a similar stripe selected 
from the 408 MHz Haslam map (which was smoothed 
to 7° and scaled to 30 GHz using /3t = -3). We de- 
tected no cross-polarization of any type between POLAR 
and the Haslam map, obtaining a mere upper limit of 
\X\ <U.OiJ,K (or a 2-cr upper limit of IbAfxR). 



III. RESULTS 

A. POLAR Power Spectra 

POLAR was a ground-based CMB polarization exper- 
iment that operated near Madison, Wisconsin [22,81,82]. 
It used a simple drift-scan strategy, with a 7° FWHM 
beam at 26-30 GHz, and simultaneously observed the 
Stokes parameters Q and ?7 in a ring of declination 
S = 43°. Because POLAR was insensitive to the unpo- 
larized CMB component, we cross-correlate their Q and 
U data with the T-data from the CO BE /DMR map [83]. 



2. Likelihood Analysis 

We complement our band-power analysis with a likeli- 
hood analysis where we assumed that B = 0. Specifically, 
we set B = Y = Z = and take each of the remaining 
power spectra (T, E, X) to be constant out to £ = 30. 
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FIG. 5. Likelihood results using the i?-polarized information 
alone (right panel, solid black line), using T information alone 
(left panel, solid black line), and using both POLAR and DMR 
T-information and marginalizing (solid rod lines). From top to 
bottom, the two horizontal lines correspond to 68% and 95% con- 
fidence limits, respectively. 

We first perform a simple 1-dimensional likelihood 
analysis for the parameter E using the POLAR data 
alone (discarding the DMR information), obtaining the 
likelihood function in excellent agreement with that pub- 
lished by [22] - see Figure 5 (right panel, solid black 
line). A similar 1-dimensional likelihood analysis for 
the parameter T using the DMR data alone produces 
T « 28yLtK, consistent with that of the DMR team [83] 
(left panel, solid black line). We then compute the likeli- 
hood function including both POLAR and DMR data in 
the 3-dimensional space spanned by (T, E, ri) and com- 
pute constraints on individual parameters or pairs by 
marginalizing as in [85]. Once again, we obtain a T- 
measurement in complete agreement with that for the 
DMR team (left panel, solid red line). 

Figure 6 shows our constraints in the (i?, )-plane af- 
ter marginalizing over T. It is seen that our constraints 
on the cross-polarization are weaker than the Schwarz 
inequality |r^| < 1, so in this sense the data has taught 
us nothing new. The likelihood function is seen to be 
highly non-Gaussian, so obtaining statistically meaning- 
ful confidence limits (which is of course uninteresting in 
our case, since the constraints are so weak) would involve 
numerically integrating the likelihood function. Since 
is expected to oscillate between positive and negative val- 
ues, using a flat (constant) ri in the likelihood analysis 
runs the risk of failing to detect a signal that is actually 
present in the data, canceling out positive and negative 
detections at different angular scales. This is not likely 
to have been a problem in our case, since ri is uniformly 
positive in our sensitivity range £ = 14 ± 5 for the con- 
cordance model. 

Figure 1 compares our results with all other polariza- 
tion constraints published to date. 




E amplii^udc [fiK] 

FIG. 6. Joint constraints on E polarization and r£ after 
marginalizing over T. From left to right, the contours show that 
the likelihood function has dropped to e~^'^, e"^" and e"*'® times 
its maximum value, which would correspond to 68%, 95% and 99% 
limits if the likelihood were Gaussian. For comparison, the concor- 
dance model predicts (E, r^)=(0. 001, 0.66) at £=14, the center of 
our window function for X (see Table 1). 

B. The Leiden Power Spectra 

1. Basic power spectra 

For the Leiden surveys, our analysis was performed 
using 10 multipole bands of width A£ = 10 for each 
of the six polarization types {T, E, B, X,Y, Z), thereby 
going out to ^ = 100. We used the Haslam map for 
the unpolarized component T, scaled and smoothed to 
match Leiden's five different frequencies. A Galactic cut 
of |5| — 25° was applied in order to match the POLAR 
observing region. We iterated the QE method once and 
chose the second prior to be a simple power law model 
consistent with the original measurement for the T, E 
and B power. The priors for X, Y and Z were set as 
zero. 

Figure 7 shows the E power spectra (top) and rg corre- 
lation coefficient (bottom) of the Leiden surveys. We find 
that all power spectra are well approximated by powers 
laws as in equation (2). The best fit normalizations A 
and slopes /? for E and B are shown in Table 2. The val- 
ues of (3 are consistent with previous analyses [4,24-29], 
showing that the slopes get redder as frequency increases. 
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FIG. 7. Power spectra of the Leiden data. Top panel shows 
the E power spectra for the five Leiden frequencies going from 408 
to 1411 MHz, while the Bottom panel shows the X cross power 
spectrum ri for two of the five Leiden frequencies (each frequency 
is represented by the same color in both plots). Intrinsic EB cor- 
relation could be present but masked by Faraday Rotation, since 
random rotations of the polarization angle would cause correlations 
to average to zero. 



Table 2 Normalization &: Spectral Index'") 
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FIG. 8. E and B likelihood values for the Leiden surveys. From 
top to bottom, and from left to right, the likelihoods are for the 
frequencies 408, 465, 610, 820 and 1411 MHz. As the survey's fre- 
quency increases, the Faraday Rotation reduces and we start to see 
a slight hint of an i?-excess. For all likelihoods, the contours corre- 
spond to 68%, 95% and 99% limits. The diagonal lines correspond 
to E = B. 



2. Is it E oris it B? 



For all Leiden surveys, the X and Y power spectra are 
found to be consistent with zero - the 2.4 GHz Parkes 
survey had a similar finding for X [29]. These are not 
surprising results: if Faraday Rotation makes the polar- 
ized and unpolarized components to be uncorrelated (see 
Figure 3), it is natural to expect that X, Y=0. However, 
at the CMB frequencies (where the effects of Faraday 
Rotation & Depolarization are unimportant) this should 
not be the case. 

To study the frequency dependence, we average the 10 
multipole bands of the Leiden power spectrum measure- 
ments together into a single band for each polarization 
type to reduce noise. From these results, we fit the av- 
erage frequency dependence (for the 25° cut data) as a 
power law as in equation (2) with slope = —1.3 and 
aB = —1.5 for E— and B— polarization, respectively. 



An interesting question about polarized foregrounds is 
how their fluctuations separate into E and B. Although 
many authors initially assumed that foregrounds would 
naturally produce equal amounts of E and B, Zaldar- 
riaga [77] showed that this need not be the case. There 
are plausible scenarios where the foreground polarization 
direction could preferentially be aligned with or perpen- 
dicular to the gradient of polarized intensity, thereby pro- 
ducing more E than B. In contrast, it is more difficult to 
contrive scenarios with more B than E, since they require 
polarizations preferentially making a 45° angle with the 
gradient. 

Early studies [25,29] have indicated that E k, B &i 
2.4 GHz in the Galactic plane. However, these analyses 
used Fourier transforms and spin-2 angular harmonic ex- 
pansions, respectively, without explicitly computing the 
window functions quantifying the leakage between E and 
B. This leakage is expected to be important both on the 
scale of the Parkes stripe thickness and on the pixel scale 
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[75,80], and would have the effect of mixing E and B 
power, reducing any E/B differences that may actually 
be present. Moreover, no study of the E/B ratio has 
ever been done on the large angular scales {i < 40), that 
are the most important for constraining reionization and 
inflationary gravitational waves. 

We therefore perform a likelihood analysis of the Lei- 
den surveys specifically focusing on this question, and 
including an exact treatment of the leakage. The likeli- 
hood analysis of the data is done with two free parame- 
ters corresponding to the overall normalization of the E 
and B power spectra, and assuming that they both have 
the same power law shape given by the slopes (3e from 
Table 2. The results are shown in Figure 8. Note that 
the E and B amplitudes are consistent with being equal 
to high accuracy at 408, 465, 610 and 820 MHz. At the 
highest frequency of 1.4 GHz, however, we see a hint of 
an E-excess at the 30% level, but this is only significant 
at a level of around 95%. This hint is intriguing, since 
it can in principle be given a natural physical interpreta- 
tion. It may be that synchrotron polarization has E > B 
at CMB frequencies, and that Faraday Rotation is hid- 
ing this underlying asymmetry at low frequencies. If the 
Faraday effect rotates each polarization angle by a for all 
practical purposes random amount, this will destroy any 
intrinsic alignment between the direction of the polariza- 
tion and the direction of the local intensity gradient and 
therefore produce equal amounts of E and B signal. 

3. Quantifying the Importance of Faraday Rotation & 
Depolarization for the CMB 

The key challenge for modeling synchrotron polariza- 
tion as a CMB foreground is to answer the following ques- 
tion: above which frequency are the effects of Faraday 
Rotation & Depolarization so small that our measure- 
ments can be safely extrapolated up to CMB frequen- 
cies? From an analysis of the Leiden surveys, Spoelstra 
[35] found an upper limit for RM of 35 rad m~^. Set- 
ting = 1 rad in equation (3), this suggests that the 
Faraday Rotation becomes irrelevant somewhere around 
2 GHz. However, considering that the determination of 
RM is poor in many parts of these surveys, this 2 GHz 
value is questionable. Moreover, because of the impor- 
tance of Depolarization which affects large scales more 
than small scales, we should expect the answer to de- 
pend on the angular scale £ considered. 

Let us now quantify this empirically. Figure 9 shows 
the synchrotron power spectra as a function of frequency 
for a sample of angular scales i. Using the fits from Ta- 
ble 2 and Equation (2) suggests that the polarization per- 
centage p = 5T^ /5TJ saturates to a constant value for 
u > 1 GHz at ^=50, > 4 GHz at £=14 and > 10 GHz 



at t=2. This suggest the following universal behavior**. 
At high frequencies, where the Faraday Rotation & De- 
polarization effects are unimportant and the polarized 
fluctuations simply constitute some constant fraction of 
the total fluctuations, we can use the same a for polarized 
and total synchrotron radiation in the CMB range. How- 
ever, moving to the left in Figure 9, one reaches a critical 
frequency v^, below which the Faraday Rotation & De- 
polarization effects suppress the polarized fluctuations. 
At this point, the power law changes asymptotes from a 
steeper (solid lines) to a shallower (dashed lines) power 
law, and the critical frequency in which this effect 
occurs change with the angular scale £. In order words, 
whether we can safely extrapolate our results up to CMB 
frequencies depends not only on the frequency but also 
in the angular scale. For instance, the contamination of 
the CMB quadrupole from Galactic synchrotron polari- 
zation can only be obtained from extrapolations of data 
at frequencies exceeding ~ 10 GHz, with i/* dropping 
towards smaller angular scales 1^^. 

All the information above is summarized in Figure 10, 
which shows contours of constant polarization percent- 
age p = 5T^ /5Tf in the two-dimensional (£, v) plane. In 
other words, this figure can be interpreted as a contour 
plot of the Depolarization. The Depolarization is seen 
to be negligible at high frequencies and on tiny scales, 
gradually increasing towards the lower left corner (to- 
wards low frequencies and on large angular scales) where 
Faraday Rotation & Depolarization effects become im- 
portant. 



**Iri the limit of high frequencies (where Faraday Rotation 
& Depolarization vanish), we expect the polarization fraction 
p to become frequency independent. It may still depend on 
angular scale £, however. If it does depend on I, there is no 
fundamental reason why it cannot exceed 100% on some an- 
gular scales (even though the polarization at a given point is 
by definition < 100%) — imagine, say, a uniform synchrotron- 
emitting plasma with small-scale variations in the magnetic 
field direction. However, Figure 7 and Table 2 show that as 
the frequency increases, the polarized power spectrum gets 
progressively redder, providing a tantalizing hint of conver- 
gence towards the same power spectrum slope as the total 
intensity compolent. It this is actually what happens in the 
high frequency limit, then the polarization fraction does in- 
deed become a simple constant. 

^^Due to the fact that we are dealing with cross-correlations, 
the results presented here should not be biased by systematic 
errors or calibration uncertainties in input data, since they 
would be uncorrelated betweeen the different surveys used. 
Spurious offsets will not cause excess noise either, since we 
removed the zero-point from each survey before calculating 
the cross-correlations. 
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FIG. 9. The B-polarized (solid) and unpolarized (dashed) 
power spectra STi of Galactic synchrotron emission are plotted as 
a function of frequency for multipoles £ = 2 (blue), £ = 14 (red) 
and £ = 50 (green) using the fits from Table 2 (corresponding to a 
25° Galactic cut data). The T curves (dashed) assume a = —2.8. 
For comparison, the POLAR upper limit of i? < TAfiK centered in 
£ ~ 14 (see Table 1) is shown in the lower right corner. Comparing 
this with the red curve implies either a low synchrotron polarization 
percentage or a steeper spectral index (lower a). 




10 100 

Multipole 1 

FIG. 10. {£, v) plane showing contours of constant polarization 
percentage. From bottom up, the curves are for the 0.01%, 0.1%, 
1%, 10% and 70%. 

This has important implications. For instance, a nice 
all-sky simulation of synchrotron polarization at CMB 
frequencies was recently performed assuming that the 
power spectra of cos26' and sin26' (where is the pola- 
rization angle) were frequency independent [29]. Our re- 
sults indicate that these two power spectra are dominated 



by Faraday Rotation & Depolarization effects, which im- 
plies that the E and B power should be mostly due to 
changes in polarization angle 6, and not to variations 
in overall intensity - this precise behavior is also seen by 
[29] . If Faraday Rotation & Depolarization effects are in- 
deed dominant, then it is not obvious that such frequency 
extrapolation of the cos 20 and sin 26 power spectra are 
valid. 

Figure 9 also shows the POLAR hmit oi E < 7.4/iK 
from Table 1 (lower right corner). Since this limit is cen- 
tered in £ ^ 14, it can be directly compared with the 
middle (red) curve. The noticeable gap between the two 
implies that we get interesting constraints from POLAR 
on foreground models. No synchrotron polarization is 
detected even though the Haslam stripe shows substan- 
tial synchrotron emission in the POLAR region, so ei- 
ther the synchrotron polarization percentage is small or 
the synchrotron emission falls even more steeply towards 
higher frequencies than the plotted curves indicate. A 
spectral index a = —2.8 (as shown in the plot) is only al- 
lowed if the polarization percentage p is lower than 10%. 
If p=20%, then a < —3.0, and almost complete polari- 
zation (about 70% is physically possible) would require 
a < —3.4, in poor agreement with theoretical and obser- 
vational indications [52-54]. In other words, our results 
suggest a rather low synchrotron polarization percentage 
at CMB frequencies^. 

IV. CONCLUSIONS 

CMB polarization and its decomposition into E and B 
modes is a topic of growing importance and interest in 
cosmology. In the era of MAP, a key issue is to estimate 
the contribution of Galactic foregrounds (more specifi- 
cally, polarized synchrotron emission) to these modes. 
We have used the POLAR experiment and radio surveys 
in order to quantify this contribution at large angular 
scales. 

Using matrix-based quadratic estimator methods, we 
cross-correlated POLAR with DMR data and obtained 
upper limits of £: < 7AfiK and |A:| < U.lfiK at 95% 
confidence. These upper limits are, unfortunately, to 
high to place intetesting constrains on reionization mod- 
els. A similar cross-correlation analysis was performed 
by replacing the DMR with the Haslam data, obtaining 
an upper hmit of \X\ < IbAfiK at 95% confidence. 



"From the CO-Bi?/DMR-Haslam cross-correlation results 
[87,88], we know that the rms Galactic signal of the syn- 
chrotron emission is lower than 7.1/iK at 53GHz. Note that 
this value is substantially lower than the one we obtain when 
extrapolating the ^ ~ 14 curve of Figure 9 to the DMR fre- 
quencies. This result indicates that a < —2.8 or that there is 
a deviation from the power law behaviour at frequences above 
a few GHz. 
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Wc also used our quadratic estimator methods to mea- 
sure the power spectra from the Leiden surveys, obtain- 
ing the following key results: 

1. The synchrotron E- and B-contributions are equal 

to within 10% from 408 to 820 MHz, with a hint of 
^-domination at higher frequencies. One interpre- 
tation is that E > B dX CMB frequencies but that 
Faraday Rotation mixes the two at low frequencies. 

2. Faraday Rotation & Depolarization effects depend 
not only on frequency but also on angular scale - 
they are important at low frequencies {v < 10 GHz) 
and on large angular scales. 

3. We must take into account Faraday Rotation & De- 
polarization effects when extrapolating radio sur- 
vey results from low to high galactic latitudes and 
from low to high frequencies. 

4. We detect no significant synchrotron TE cross cor- 
relation coefficient (|r| < 0.2), but Faraday Rota- 
tion could have hidden a substantial correlation de- 
tectable at CMB frequencies. 

5. Combining the POLAR and radio frequency re- 
sults, and the fact that the i?-polarization of the 
abundant Haslam signal in the POLAR region is 
not detected at 30 GHz, suggests that the syn- 
chrotron polarization percentage at CMB frequen- 
cies is rather low. 

Experiments such as MAP and Plauc;k will shed signif- 
icant new light on synchrotron polarization and allow 
better quantification of its impact both on these experi- 
ments and on ground-based CMB observations. 
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